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Abstract 
>v 

We present the class of regular homogeneous T-models with vacuum dark fluid, associated 
■ with a variable cosmological term. The vacuum fluid is defined by the symmetry of its stress- 

\^ \ energy tensor, i.e., its invariance under Lorentz boosts in a distinguished spatial direction (pj = 

7— I ' — p), which makes this fluid essentially anisotropic and allows its density to evolve. Typical 

features of homogeneous regular T-models are: the existence of a Killing horizon; beginning 
of the cosmological evolution from a null bang at the horizon; the existence of a regular static 
pre-bang region visible to cosmological observers; creation of matter from anisotropic vacuum, 
^jpj 1 accompanied by very rapid isotropization. We study in detail the spherically symmetric regular 

l — 1 . T-models on the basis of a general exact solution for a mixture of the vacuum fluid and dustlike 

I \ matter and apply it to give numerical estimates for a particular model which illustrates the 

J> ■ ability of cosmological T-models to satisfy the observational constraints. 

00 : 

\o . 

! 1. Introduction 

Astronomical data provide convincing evidence that our Universe is dominated at above 70 per cent 
of its density by a dark energy responsible for its accelerated expansion due to negative pressure, 
p = wp with w < —1/3 [1-5]. Observations constrain the equation of state of dark energy within 
\ —1.45 < w < —0.74 [4], the best fit w = — 1 [6-11] corresponds to a cosmological constant A related 
^ to the vacuum density p vac by A = 8nGp vac (for a detailed discussion see [4,11,12]). 

The Einstein cosmological term Ag^ u = S7iGp vac g^ u is plagued with the cosmological constant 
problem. Quantum field theory estimates A by the Planck scale [13], so that the resulting zero-point 
density is incompatible with all observational data [14], and a long-term problem was how to zero it 
out [13]. A key problem now is why the cosmological constant should be at the scale demanded by 
observations [12]. Another dynamical aspect is that the inflationary paradigm needs a large value of 
p vac at the earliest stages to provide a reason for rapid initial expansion [15], typically at the GUT 
scale [16] (for a review see [17]), observations testify that it is small at the present epoch, whereas 
the Einstein equations require p vac = const. 

A lot of theories and models have been developed to describe dynamically evolving dark energy 
which today behaves like a cosmological constant (for a recent comprehensive review see [12]). 

On the other hand, the Einstein equations admit a class of solutions with source terms associated 
with vacuum denned by the symmetry of its stress-energy tensor, including vacuum with reduced 
symmetry which allows the vacuum density to be evolving and clustering [18,19]. 

The Einstein cosmological term Ag^ is associated with the maximally symmetric vacuum stress- 
energy tensor (all eigenvalues equal) T^ v = p^g^ [20] representing de Sitter vacuum with p = — p vac 
and p vac = const. The symmetry of a vacuum stress-energy tensor can be reduced to the case 
when only one or two of the spacelike eigenvalues of T^ v coincide with its timelike eigenvalue, 
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Tj = T\ [18,21,22], so that the vacuum equation of state holds in the j-th direction(s), pj = —p, 
and thus a vacuum with reduced symmetry must be evidently anisotropic. Evolution of vacuum 
density and pressures is governed by the conservation equation T^ v v = 0. 

The class of GR solutions specified by pj = —p, which describe time-dependent and spatially 
inhomogeneous vacuum energy, presents, in a model-independent way, anisotropic vacuum fluid 
with continuous density and pressures. It can provide a unified description, based on a space- 
time symmetry and introduced in [22], of dark ingredients in the Universe by a vacuum dark fluid 
representing both distributed vacuum dark energy (by an evolving and inhomogeneous cosmological 
term [18,23]) and compact self-gravitating objects with interior de Sitter vacuum [19,24-26]. 

The variable cosmological term A^ v = 87rGT At ^( vac ) has been introduced in the spherically sym- 
metric case [18] when a vacuum T^^az) is invariant under radial Lorentz boosts [21], while the full 
symmetry is restored asymptotically at the centre and at infinity [19,21,24,27]. 

It provides a description of a smooth evolution of the vacuum from Ag^ u to \g^ u with A < A . In 
the spherically symmetric case, T^^) generates globally regular space-times with a de Sitter center 
[19]. Depending on the parameters and on the choice of a reference frame, they describe localized 
objects with de Sitter vacuum trapped inside (nonsingular black and white holes [18,21,28,29], 
self-gravitating vacuum structures without event horizons [24] called G-lumps [19,22]), and regular 
cosmological models with variable vacuum density and pressures considered and classified in our 
previous paper [23]. 

In the cosmological context, the spatial direction x- 7 in which the vacuum symmetry of A^ IV = 
87rGT At j y ( vac ) is preserved, is necessarily distinguished by the expansion anisotropy. In cosmologies 
guided by A^ u , an anisotropic stage is generic [23,28]. A further early evolution is the most 
transparent in homogeneous T-models, and this is the question addressed in the present paper. 

T-models are generically related to the T-regions found originally by Novikov as essentially non- 
static regions beyond the horizon in Schwarzschild space-time [30] and specified in the spherically 
symmetric case, ds 2 = e v ^ q '^dt 2 — e^ q '^dq 2 — r 2 (q,t)dQ 2 , by the invariant quantity A = g %k r^r^ = 
e~ u (dr/dt) 2 — e~ fJ '(dr /dq) 2 [31,32]. When A > 0, the normal Ni = r t i to the surface r = const 
is timelike, so that the surface r = const is spacelike, and orbits of constant r and hence static 
observers do not exist in principle. T-regions are further specified as T_ when f < and T + when 
r > (with A > 0, f cannot vanish, and this condition is invariant) [32] 

Formally, spherically symmetric T-models were introduced by Ruban [35] as modification of vac- 
uum T-regions to the case of dustlike matter considered in a synchronous comoving reference frame 
with the metric representing an inhomogeneous generalization, ds 2 = dt 2 — e a<yt ' q ^dq 2 — R 2 {t)dVt 2 , of 
the anisotropic cosmological models of semi-closed type [36] with hyper cylindrical spatial sections 
V3 = §2 ® Actually Ruban deeply revised and generalized (to the case of nonzero A) Datt's 
solution [37] qualified by the author as being "of little physical significance" [38]. The particular 
case when a(t,q) = a(t) is known as the Kantowski-Sachs model [39]. A remarkable feature of 
T-models discovered by Ruban is a gravitational mass defect maximally possible in GR, equal to 
the total rest mass of the matter which appears to be entirely gravitationally bound [35]. 

T-models were further extended to the case of a perfect fluid with non-zero pressure [40-44], and 
the class of spherically symmetric T-models with isotropic perfect fluid was completed by Ruban in 
1983 [45]. 

Amazingly, Ruban excluded from his consideration the case Tq = Tf (p r = — p) as being "of 
no interest for cosmology and relativistic astrophysics" [46] . On the other hand, the requirement of 

1 Thcre are exceptions from this rule, represented by T-regions with r(t) having a minimum. This happens, e.g., in 
regular black-hole solutions of brane gravity [33] and in some recent solutions with nonminimally coupled Yang-Mills 
fields [34]. 
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regularity leads in this case to globally regular space-times with an obligatory de Sitter centre [19] 
which represent distributed and clustered vacuum dark energy and describe both black (white) 
nonsingular holes and regular cosmological models dominated by anisotropic vacuum fluid, including 
regular modifications of Kantowski-Sachs models [23]. In these models, the cosmological evolution 
starts from horizons with a highly anisotropic null bang, but information on the pre-bang history 
in the de Sitter region in the remote past is available to Kantowski-Sachs observers [23,47]. 

The latter circumstance provides a solution to one more long-standing problem, that of the 
initial cosmological singularity. 

Homogeneous regular cosmological T-models with anisotropic vacuum fluid can be seen as T- 
regions of space-times with 3-parameter isometry groups acting transitively on spatial 2-surfaces, i.e., 
those with spherical, planar and pseudospherical symmetries. Among them spherically symmetric 
T-models related to Kantowski-Sachs (KS) type regular models are of greatest interest. 

The aim of this paper is to introduce and to study spherically symmetric homogeneous regular 
cosmological T-models with anisotropic vacuum fluid and dustlike matter, bearing in mind that 
the material content of the modern Universe is almost entirely represented by dark energy and 
pressureless matter. 

The paper is organized as follows. In Sec. 2 we consider T-models of KS type filled with dust and 
an anisotropic vacuum fluid described by A Mi/ . For non-interacting dust and vacuum, we present 
a general solution without specifying the particular density profile of A^ u . In Sec. 3 we argue that 
regular models on the basis of the above solutions must have a static, purely vacuum, asymptotically 
de Sitter region in the remote (pre-bang) past, and briefly describe the properties of such regular 
models. In Sec. 4 we give a detailed description of a particular model illustrating the ability of 
regular T-models to satisfy the observational constraints. In Sec. 5 we summarize and discuss the 
results. 



2. General solution 

We begin with the KS metric in the general forrr|§ 

fa* = e 2 ldt 2 _ e 2* dx 2 _ e 2P dQ 2 = ^2 _ ^^2 _ ^^2 (1) 

where a, (3, 7 are functions of the time coordinate t and dfl 2 = dO 2 + sin 2 Odtp 2 . Then the nonzero 
Ricci tensor components are 

R\ = e" 27 [q + 2/3 + a 2 + 2/3 2 - j(a + 2/3)] ; (2) 
R% = e- 27 [a + d(d + 2/3- 7 )]; (3) 
R 6 Q = R%= e~ 2fS + e" 27 [/3 + $(a + 2$- 7)] (4) 

(the dot stands for d/dt). Two Einstein tensor components to be used are 

G\ — — e~ 2f3 - e-^0 2 + 2aj3), (5) 
G x x = - e" 2/3 - e" 27 (2/3 + 3{3 2 - 2/3 7 ). (6) 

The Einstein equations read 

g; = r;- = -sttGt;. (7) 



2 Our conventions are: the metric signature (H ) ; the curvature tensor R" ^ pi , = d u TZ — . . . , R^ — R" ' , 

so that the Ricci scalar R > for de Sitter space-time and the matter-dominated cosmological epoch; the system of 
units c = % = 1 . 
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For our purposes, it is helpful to choose the time coordinate t in a way similar to choosing the 
curvature coordinates in static spherical symmetry, namely, to put t = r = e 13 , so that the metric 
is rewritten as 

ds 2 = e 2 ^dr 2 - e 2a ^dx 2 - r 2 dQ 2 . (8) 
Then the (j) component of and the difference (°) — (*) read 

I^e-^ + l)]' = 8nGT*, (9) 



e~ 2 \a' + i) = 8ttG(T; - T£), (10) 



-7 + 4vrG J r e 2 \T r r - T x x )dr. (12) 



where the prime denotes d/dr. In full analogy with static spherical symmetry, Eqs. ([9]) and ([10 
may be rewritten in an integral form: 

e _ 27 = _ 1 + 8vrG [ T * r 2 d ^ (1V 
a ■ 

If T£ is known as a function of r , then Eq. (11 II) is a solution to (Q that yields the lapse function 
e 27(r) ^ anc j £ ne other metric function a is found from ffl2l) if, in addition, the difference T r r — T£ is 
known as a function of r . It is really so in many important cases. 

Now, let us take the stress-energy tensor (SET) of matter as a sum of SETs of dust and vacuum. 
Dust, with the density Pd(r), is at rest in a T-model with the metric (PQ), while for a spherically 
symmetric vacuum fluid we have T^( vac ) = diag(pt,, p„, — p v ±, —p v ±). The full SET reads 

= di&g(p d + p v , p v , -p v ±, -p v ±). (13) 

Evidently, if we know the vacuum density profile p v (r) , Eq. (1111) gives us e 27 ^: 

-27 = _! + M{r) :=A7rG [ p v r 2 dr. (14) 
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The function A^(r) is similar to the mass function conventionally introduced in the description of 
static, spherically symmetric configurations. 

Note that this result is obtained without using the conservation law for dust and vacuum and 
is even valid when they are interacting. 

It remains to find the dust density pd and the second metric function a(r) . This can be done 
with the aid of the conservation law component V M T^ = which in our case reads 

p' d + p d (a' + -) +p' v + ^(p v +p v ±) = 0, (15) 

and Eq. ( TlOl) where T r r — T£ = pd- At stages of the cosmo logical evolution when there is no interaction 
between dust and vacuum, we can write the conservation law separately for each of them. Then, 
for dust we have 

p' d + Pd (at' + = Pd = Poe~ a /r 2 , p = const, (16) 

while the corresponding relation for vacuum 

f/ v + -(Pv+Pv±) = (17) 
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expresses p v ± in terms of p v . Substituting ([TBI) into (jTUI) . we can rewrite it as follows: 




(18) 



which after integration gives 



AttGhq e 



/ 




(19) 



e 7 being given by ( fl4l) . This completes the solution. The physical time r is obtained by integration, 



and all unknowns are expressed as functions of r in a parametric form. 

Thus we have solved by quadratures the Einstein equations for KS cosmology with a non- 
interacting mixture of dust and -vacuum. The solution contains one essential integration con- 
stant /io and one arbitrary function, the vacuum density profile p v {r) (or equivalently the mass 
function Ai(r)). This is a new solution for the general case p v ± ^ —p v - In the particular case of 
the conventional cosmological constant A = 8irGp v = —8nGp v ± = const, our solution agrees with 
Ruban's homogeneous solution [35]. 

For interacting dust and vacuum, a further solution depends on the particular form of the 
interaction. Though, if we know (or have a reason to prescribe) the function pd(r), then Eq. ([T2~]) 
immediately yields a as a quadrature. 

3. Nonsingular models and the static core in the pre-bang past 

Let us now discuss the possibility of obtaining globally regular models of the Universe on the basis 
of the above KS metrics. The following three points will be argued: 

(i) In regular expanding Universe models, the cosmological evolution should begin from a Killing 



(ii) The region beyond the horizon should be purely vacuum, i.e., dust is absent there. 

(iii) Dust appears in the cosmological region due to interaction with the vacuum fluid. 

To prove item (i), we will use, in the metric (pj, a new time variable specified by the coordinate 
condition a + 7 = 0, so that (TjJ takes the form ds 2 = a~ 2 {t)dt 2 — a 2 {t)dx 2 — r 2 (t)dQ 2 . Then, instead 
of ( [TOl) . we obtain the equation 



Assuming p^ > 0, we immediately conclude that r < 0. Consequently, if there is expansion (f > 0) 
at some instant ti, we can assert that r(t) reached zero at some finite earlier instant t s < t±. 
Meanwhile, r = is a curvature singularity of the metric ([TJ. 

The only way out is to assume that the cosmological evolution started later than the would-be 
singularity, which is only possible if there was a Killing horizon at some th > t s , where a(th) = 
0. Then there is a hope to find a static regular region beyond such a horizon. This reasonable 
assumption is additionally justified by the consideration below concerning (ii) and by the fact that 




(20) 



horizon. 



2a 2 - = -8nGp d . 



(21) 



r 
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all related regular vacuum models without dust [23] do have horizons and regular R-regions with a 
de Sitter core in their remote past. 

As to item (ii), there is a serious physical argument. If we imagine that cosmological dust is 
continued to the static region where the cosmological time t becomes a radial coordinate and vice 
versa, then the density pd = T*, d , is converted there to radial pressure, while the pressure —T^ = 
is converted to zero density, and we obtain quite an unplausible kind of matter with zero density 
but nonzero pressure. 

Moreover, it can be strictly shown that pd = at the horizon. Indeed, at t = th, according to 
(ED), Pd = due to a = 0, unless r — > oo as t — > th- But the latter opportunity is unacceptable 
for a horizon as a regular surface where r 2 (t) should be an analytic functional Thus pd = at the 
horizon, and it is natural to conclude that dust did not exist to the past of it (item (ii)). 

Returning to the cosmological region, we notice that non-interacting dust is incompatible with 
a horizon: it would then have an infinite density according to ( Fl6|) . pd ~ [a(t)r 2 (t)] _1 . Since pd = 
at the horizon, pd > at later times can only emerge from the vacuum fluid due to interaction - 
and this is our item (iii). 

Let us now consider the region beyond t = th filled with vacuum fluid and suppose that there is 
a regular centre r = 0. Such space-times were discussed in detail in Ref. [23]; they can be described 
by metrics of the form ([8]) with e -27 = e 2a = —A(r) (so that r is a spatial coordinate while a; is a 
temporal one), i.e., 

dv 2 

ds 2 = A{r)dtl ~A(r) rW ' (22) 
with A(r) = — e~ 27 given by Eq. ffl4]F l. Taking at the centre A(0) = 1 , we can write 

A( r ) = 1 - y^LL^ M(r) = AttG f p v r 2 dr } (23) 
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so that M. is the conventional mass function. 

The function A(r) may have different number and orders of zeros (depending on the particular 
density profile p v (r)), corresponding to Killing horizons in space-time [23]. 

In any space-time with the metric (1221) . the geodesies equations have the following first integral: 



—J + kA(r) + —A(r) = E 2 , (24) 



where E and L are the constants of motion associated with the particle energy and angular momen- 
tum, while the constant k takes the value k = 1 for timelike geodesies and k = for null geodesies; 
the afline parameter r has the meaning of proper time along a geodesic in the case k — 1 . At a 
horizon r = r^, A(r) vanishes, and in case E / one has dr/dr ^ for all geodesies, whence it 
follows that |r| < oo, irrespective of the order of the horizon. For KS comoving observers whose 
trajectories are time lines in the T-region, A = 1, L = 0, E 2 = 0. Therefore even "slow" observers, 
whose world lines coincide with these time lines, receive information from their infinitely remote 
past brought by particles following the geodesies (124")) . So, KS observers, even in the case of a double 
horizon and a null bang in their infinitely remote past, can receive pre-bang information brought by 
particles and photons following other geodesies and crossing the horizon at their finite proper time. 



3 The coordinate defined by the condition a + 7 — varies near a horizon (up to a nonzero constant factor) like 
manifestly well-behaved null coordinates of Kruskal type used for analytic continuation of the metric [54,55]. This 
implies the analyticity requirement for r 2 (t) . 

4 Note that the time coordinate of the static region has nothing to do with the cosmological time in the metric 
(TT]) : it is rather related to the spatial coordinate x . 



7 



Let us here discuss the simplest generic structure with a single simple (first-order) horizon, so 
that the global space-time structure is like that of de Sitter space. The horizon radius r = 
corresponds to = 2A / f(r/ l ) and can be found if p v {r) is specified. 

The static R-region between r = and r = r% represents a regular "core" in the remote past of 
an expanding KS universe in which the cosmological evolution starts with a Null Bang (i.e., contains 
a horizon instead of a singularity). 

It is important that Eq. (1141) with Ai(r) given in (123]) is now valid for all values of r, in the whole 
space-time, whereas for the other metric coefficient e 2a there are different expressions: e 2a = —A{r) 
for r <rh, (compare (E]) and (|22l); furthermore, Eq. ( fl9l) holds for the cosmological epoch without 
interaction, and there is a yet unknown expression for e~ 2a in the interaction epoch. 

Since the horizon is crossed at a purely vacuum stage, a nonzero matter density pd appears due 
to interaction between dust and vacuum in a certain period after the Null Bang. This perfectly 
conforms to the well-known prediction of QFT in curved space-time that nonstationary space-times, 
especially anisotropic ones, create particles [36,48,49], and a decaying vacuum density accompanied 
by growing pd is a phenomenological description of this process. 

In what follows we will give some estimates for a specific choice of the vacuum decay law. 

3/3 

4. T-model with vacuum density p v = p c e~ r ' r * + p\ 

Let us choose the profile p v (r) in the form 

Pv = p c e~ r3/r * + pa, Px = r^, p c , r*, A = const > 0, (25) 

The first term results from a simple semiclassical model for vacuum polarization in spherically 
symmetric gravitational fields [19,21,24]. The second one refers to a background cosmological 
term which is at the scale suggested by observations. An additional indirect justification of the 
ansatz ( 1251) comes from the minisuperspace model of quantum cosmology. In a certain gauge, the 
cosmological constant is quantized as an eigenvalue of the appropriate potential in the Wheeler- 
DeWitt operator [50]. The form of the potential corresponds to a non-zero quantum value of the 
cosmological constant at each finite value of the time-dependent scale factor. 

We thus have three distinct length scales: r c = y^3/(8nGp c ) , a de Sitter radius corresponding 
to the central density p c , the scale r* ~ (r 2 ?-^) 1 / 3 [21,51] characterizing the vacuum decay rate 
in the spherically symmetric gravitational field (r g is the Schwarzschild radius related to the total 
gravitational mass of the decaying vacuum), and the scale r A related to the background cosmological 
constant A, assumed to be of the present Hubble order of magnitude, A ~ 10 -56 cm -2 , r\ ~ 10 28 
cm. 

Let us assume that the central density p c is of the GUT scale, 
p c w (10 15 GeV) 4 w 2.2 • 10 75 g/cm 3 r c w 0.8 ■ 10" 25 cm, (26) 

and r* rs 4.6 • 10~ 8 cm for the scale r* given in [21]. 

Horizon radius. Since r* ^> r c , evidently, the density p v is nearly constant and the metric is 
approximately de Sitter at radii r < r c , and, with good accuracy, we have a cosmological horizon at 
r Rj r c . The same is true for any profile p v (r) with a slow enough decay rate from p c corresponding 
to a regular centre. As a result, in all such cases, including ([25]) . 



Pv{r h ) Ri p c . 



(27) 



s 



More precisely, with (1251) and (|26|) . we have = r c [l + 0(6 3 )] where 6 = r*/r c P3 5.6 • 10 17 . 

Vacuum decay time. Next, let us estimate the vacuum decay time, taking r = at crossing the 
horizon and seeking t/, the instant when the first, exponential term in fl25l) becomes equal to the 
second term, p\ . So, 

T f =[ f e 7(r) cir, (28) 



where r/ is found from the condition 

p c e- r /^=p A r/«6.3r*. (29) 
The integrand e 7 is expressed from the relation due to ( 1231) 

e" 27 = -1 + -(1 - e"* 3 ) + -r 2 x 2 , x := -. (30) 
x 3 r* 

The integral (1281) is taken over the interval (1/6, 6.3) . It is easily seen that, on this interval, the 
term containing A is more than 60 orders of magnitude smaller than the others and can be omitted. 
Near the horizon, where x <§C 1 , the integral may be calculated analytically: 



n=7- / [r,-rj« J' 



dx 



2 x 2 -l 

where r\ is some small enough value of r. For r\ = 0.1 r* we have (with a relative error smaller 
than 10 -3 ) 7~i ~ 39.3 r c . In finding the remaining integral, we may leave in Eq. (l30l) . with high 
precision, only the term with 6 2 and get numerically 

f 63 J^dx 
r 2 =r c / 3 ^12.3r c , 

J o.i v 1 - e x 
which finally gives 

T/ = n + r 2 ss 51, 6 r c ss 0.4 ■ 10" 23 cm, (31) 

or, in seconds, ry P3 1.3 • 10~ 34 s. This value is rather close to the GUT parameters, while by this 
time the scale factor r has already inflated from r = P3 0.8 • 10~ 25 cm to rj ^ 6.36 r c f=3 2.8 • 10 -7 
cm. The other scale factor, a(r) , has inflated from zero at the horizon to some finite value. 

The inflation of r corresponds to approximately 43 e-foldings which is not regarded sufficient 
in conventional inflationary cosmology. However, in our class of models inflation as such is not 
necessary because, due to the existence and observability of a static core in the remote past, all 
parts of our model Universe are causally connected. 

Solution in the regime of constant A. In the period r > t/ we can neglect the variable part of 
the vacuum SET, remaining with the usual cosmological term, p v = —p v ± = \/(8ttG) > 0. Then, 
according to ([HP , the lapse function is given by 

^ 2M 1, „ 

e" 27 = -1 + + -Ar 2 , 32 

r 3 

where M is an integration constant expressing the total contribution of the decaying vacuum 
component to the mass function (this corresponds to dropping e~ x in ( 130]) . and we have then 
2M = 6 2 r* P3 1.4 • 10 28 cm). 
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Integration in Eq. (|T9|) then leads to a very cumbersome expression in terms of elliptic functions, 
and we will not present it here. We notice that, in the absence of dust, the integral in (Tl9|) reduces 
to an arbitrary constant, and after simple rescaling of t (or directly from Eq. (flOl) ) we obtain the 
conventional Schwarzschild-de Sitter metric in its T region, i.e., a = —7 with 7 given by 



Matter-dominated epoch. In our model, there is a large period when 

r/ < r < M ~ 10 28 cm, (33) 

in which ^> p v and we can write e 27 ~ 2M/r. This corresponds to the epoch when all matter 
has been already created but the influence of the cosmological constant A is still negligible. In this 
period, the integral in ffl9l) is easily calculated, and both r and a may be expressed in terms of the 
cosmic time r: 



8nGfi ( Ci 



<r)= ea = ^r( r + 4)> r=[3^M/2(r-r )\ , (34) 



2/3 



where C\ and r are integration constants. This very simple expression gives a reasonable approxi- 
mation to the corresponding exact dust solution [45]. 

The values of C\ and r are determined by the details of the interaction period r < 77. If we 
admit for certainty that the matter-dominated epoch begins at r = 77 and, as before, put r = 
at the horizon, we must require r{rf) = 77 , which leads to 17 — tq ~ 0.8 • 10~ 24 cm, and using ( 1281) 
we obtain tq »s 3.2 • 10~ 24 cm. 

The value of c\ is severely restricted by the conditions a > and a' > at r > 77: 

-rj /2 < ci < 2rJ /2 . (35) 



Isotropization. The condition that a KS universe expands isotropically is that the directional 
Hubble parameters, defined as [52] 

H a := e~ 7 d = da/dr, and H r := e~ 7 /3 = df3/dr = (l/r)(dr/dr), (36) 

are equal. Entire isotropy would mean that R* = R e e , which, according to (J3]) and (@J, would lead 
(in the proper time gauge t = r) to 

a + 3d 2 = a + 3d 2 + 1/r 2 . 

This shows that any KS cosmology is necessarily anisotropic for topological reasons and can only 
isotropize asymptotically at large r. 

The degree of anisotropy may be characterized by the ratio A = a/H where 

a = (H a - H r )/VS, H=(H a + 2H r )/3 (37) 

are the shear scalar and the mean Hubble parameter, respectively. Since H a and H r can in principle 
be equal, the observable parameter A can be zero despite the topological anisotropy. 

Observations of the cosmic microwave background show that our Universe is highly isotropic [53] : 
at the recombination epoch, at which the electromagnetic radiation has decoupled from matter (at 
redshifts z ~ 1000), the ratio a/H could not exceed 10~ 6 , whence it follows that at present it is at 
most of order 10 -9 [52]. This, in principle, strongly constrains the parameters of any cosmological 
model based on the KS metric. We shall see, however, that our model is automatically sufficiently 
isotropic. 
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In terms of the quantities r and a(r) , the anisotropy parameter reads 

*=*W£- (38) 

For the matter-dominated epoch we obtain from (I34p 

A = ~2r^+ Cl - (39) 

From the restriction (1351) it is evident that our model Universe expands almost isotropically very 
soon after the end of the interaction epoch, and at r ^> rf we have simply 

A < {T f lrfl\ (40) 

For instance, in the recombination epoch (r ~ 10 25 cm), still belonging to the matter-dominated 
period, we obtain A < 10~ 48 , a degree of isotropy far exceeding the observational constraint. 

Matter density after creation. It is reasonable to suppose that, by the time t — Tf, when 
the rapidly decaying vacuum component almost completely vanishes, the whole amount of matter 
particles had already been created, and the further evolution proceeded with each of the conservation 
laws f|T7|) and ffl~6]) for dust and vacuum being valid separately, hence we have the expression f TT9]) 
for a(r) = e a (r) . 

Let us obtain an order-of-magnitude estimate of matter density at r = Tf , using the conservation 
law (fl6|) . assuming that at present, when we can put r = a = r$ ps 10 28 cm the dust density 
corresponds to the total (visible plus dark) matter density of the Universe, i.e., 

Pd ~ m /rl ps 2 • 10 -30 g/cm 3 , m = const ps 2 • 10 54 g. (41) 

Neglecting the anisotropy (which, as we saw, is really negligible soon after t — Tf), i.e., assuming 
a(r) ^ r at all r > Tf, we obtain 

p d (T f ) ~ 0.8 ■ 10 74 g/cm 3 ~ Pc . (42) 

This seems to be a reasonable figure for the state of the Universe right after the vacuum decay and 
matter creation. 

5. Conclusion 

A general solution for homogeneous T-models has been found for a mixture of vacuum dark fluid 
and dustlike matter. It presents the class of T-models specified by the density profile of the vacuum 
fluid p v . The solution contains one arbitrary integration constant related to the dust density. 
We have proved the existence of regular homogeneous T-models with the following structure: 

(i) A regular static R-region in the pre-bang past, which is asymptotically de Sitter as r —>■ 0. 

(ii) Killing horizon(s) separating static R-region(s) from KS T-region(s). 

(iii) A null bang from a Killing horizon, followed by matter creation from an anisotropic vacuum, 
accompanied by rapid isotropization. 



(iv) Further evolution in the regime without interaction between matter and vacuum. 
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Let us emphasize that the rapid isotropization was not imposed by hand. For our regular 
T-models rapid isotropization is generic, since the cosmological evolution starts from the Killing 
horizon, in this case a > 0, a > which immediately leads to tight constraints on the anisotropy 
factor. 

The most interesting feature is that information about pre-bang history is available, and hence 
observable for the comoving observers in T-regions and can influence dynamics. 

Let us also emphasize that regular homogeneous T-models with anisotropic vacuum dark fluid 
do not not need an inflationary stage since there are no causally disconnected spatial regions. Still, 
in this class of models, there is a certain amount of inflation, which happens simultaneously with 
isotropization, right after the Null Bang: the scale factor a(r) grows from zero to a finite value 
while the other scale factor r(r) rapidly grows from one finite value to another being driven by the 
(currently large but decaying) variable A-term. 
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